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1. SCOPE OF THE PAPER 

Binary stars are home to so many different physical processes that an exhaustive review 
of them would fill a thick textbook. It is therefore necessary to delineate right away the 
scope of the present text, which focuses on the methods to detect binaries. Besides the 
traditional methods involving radial-velocity or photometric monitoring, the discussion 
covers as well less known methods involving astrometry or maser (non-)detection. For 
a more traditional review focused on the evolutionary aspects (including a detailed 
description of the evolution of orbital elements), we refer the reader to the recent book 
by Eggleton yj], or to the older Saas-Fee course 

Because of the authors' own biases, this review on binary detection methods must be 
understood in the general scientific context of systems involving low- and intermediate- 
mass ('L&IM') stars, i.e., stars which end their lifes as white dwarfs (WDs). For this 
reason, some of the detection methods discussed in the present review only apply to 
L&IM binary systems; conversely, methods that are specific to massive binaries (Wolf- 
Rayet binaries, High-mass X-ray binaries, binaries involving black holes or pulsars) will 
not be addressed here. 

This review is organized as follows. The specificities of L&IM stars are briefly 
discussed in Sect. [21 Basic concepts about binaries (orbital elements and Roche 
lobe) are summarized in Sect. [31 Possible ways to detect binary systems (some- 
times specific to L&IM binaries) are reviewed in Sects. [H to [71 covering spec- 
troscopic methods (Sect. [H), photometric methods (Sect. [51), astrometric methods 
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(Sect. [61) and miscellaneous other methods (Sect. |7]). We stress the difficulties asso- 
ciated with these methods. The reader is referred to a recent lecture by Halbwachs 

(|http : / /www, astro . lu . se/ ELSA/pages/PublicDocuments /Halbwachs .p df ) for a 

good introduction on the technical aspects of the main methods, and to the book Ob- 
serving and Measuring Visual Double Stars [3] for methods relating to visual binaries, 
which are not addressed here. A summary of what may be known about masses for 
the different kinds of binaries is presented in Sect. [H and finally Sect. |9] provides an 
extensive list of internet resources (mostly catalogues/databases of orbits and individual 
measurements) for the study of binary stars. 



2. THE FAMILIES OF L&IM BINARIES 

What makes L&IM stars so interesting is that in the course of their evolution they go 
through the asymptotic giant branch (AGB) phase of evolution. AGB stars have two 
very important properties: they are home to a rich internal nucleosynthesis and exhibit 
strong mass loss. Besides controlling the AGB evolutionary timescale, the strong mass 
loss has important side effects when it occurs in a binary system, like the development of 
symbiotic activity, of associated X-ray emission, of maser emission (possibly suppressed 
by the perturbation induced by the companion), etc... In many 'after-AGB' system^ 
the mass transfer from an AGB star has left its mark on the companion, enhancing its 
abundances with the products of the AGB nucleosynthesis, most remarkably C, F, and 
elements heavier than iron produced by the s-process of nucleosynthesis [see la for a 
recent review]. An exemplary case of after-AGB systems are barium stars: G-K type 
giants remarkable for their overabundances of Ba [7] . Related families include the Abell- 
35 subclass of planetary nebulae [8J, barium dwarfs [including the so-called WIRRing 
stars; @], subgiant and giant CH stars 113, [O], extrinsic S stars, as opposed to intrinsic 
S stars, which exhibit spectral lines of the element technetium, a product of s-process 
nucleosynthesis which has no stable isotopes [12] and d'-type yellow symbiotics (^^. 
But not all of the after-AGBs need to be s-process rich. The post-AGB binaries are 
an interesting case, as they are all by definition after-AGBs: some of them do exhibit 



s-process enhancement while others do not Ill4l ll5h. Neither do symbiotic stars (SyS) 
involving M giants and massive WD companions {Mwd > 0.5 Mq) exhibit s-process 
enhancements [16]. Finally, some of the cataclysmic variables with massive WDs should 
also belong to the after-AGB family. An extensive list of families of binary stars with 
WD companions is presented in the dedicated reviews of Jorissen [|4t] and Parthasarathy 
et al. fl7]. How the different families mentioned above fit in a coherent evolutionary 
scheme is sketched in Fig.[T] 

While Fig. [D provides a classification of L&IM binaries in terms of a temporal 
sequence, one may also try to order them in terms of their physical properties, like 



The term after-AGB binaries is used here to refer to binary systems where at least one component has 
gone through the AGB. After-AGB systems should not be confused with the more restricted class of post- 
AGB systems, denoting the short transition phase between AGB and planetary nebula stages of (single or 
binary) stellar evolution. 
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FIGURE 1. The evolution of a system consisting initially of two L&IM main-sequence stars. The left 
column corresponds to the normal evolutionary sequence of single stars, while the right column represents 
the various classes of stars with chemical peculiarities specifically produced by mass transfer across the 
binary system. Hatched circles denote stars with atmospheres enriched in carbon or heavy elements. (From 
il) 



metallicity, spectral type and symbiotic activity. Symbiotic activity is expected when 
a compact star is heated by matter falling onto it from a mass-losing (often giant) 
companion. As it will be discussed in more details later on, the hallmark of SyS is a 
hot spectral continuum superimposed on cool spectral features, but symbiotic activity 
also means outbursts. X-ray emission, high-excitation emission lines, and nebular lines. 

Fig.[2lis an attempt to classify L&IM binaries in the plane metallicity - spectral-type, 
and to correlate this location with the presence or absence of symbiotic activity [for 
reviews, see[5l [l6l [l8ll. 

Metallicity, plotted along the vertical axis, has a strong impact on (i) the spectral 
appearance, which controls taxonomy (CH giants for instance - box 1 - are the low- 
metallicity analogs of the barium stars - box 6); (ii) the efficiency of heavy-element 
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FIGURE 2. The various kinds of symbiotic stars (SyS) and of peculiar red giants, in a plane spectral- 
type vs. metallicity. Within each box, the first line (slanted font) lists the stellar family, the second line 
(calligraphic font) indicates whether or not stars from that family are enriched in s-process elements 
(either from internal nucleosynthesis - "Intrinsic Ba" - or from mass transfer across a binary system - 
"Ba" standing for Extrinsic Ba). The last lines (regular font) provide the binary properties of that family: 
binary (short or long periods) or non-binary, SyS or non-SyS. Each box has been assigned a number, for 
easy reference in the text. The horizontal line with an arrow in the middle of the figure is to indicate that 
the thermally-pulsing AGB phase (where the s-process nucleosynthesis takes place) involves different 
spectral types at low- and near-solar metallicity. (From [jst]) 



synthesis, being more efficient at low metallicities JT^, and (iii) the location of evolu- 
tionary tracks in the Hertzsprung-Russell diagram (hence the correspondence between 
spectral type and evolutionary status, like the onset of thermally-pulsing AGB, where 
the s-process operates, will depend on metallicity). Fig. [21 therefore considers three dif- 
ferent metallicity ranges: (i) [Fe/H]< —1, corresponding to the halo population; (ii) 
— 1 <[Fe/H]< 0, or disk metallicity; (iii) [Fe/H]> 0, solar and super-solar metallicities 
found in the young thin disk. 

The horizontal axis in Fig. [21 displays spectral type. At a given metallicity, spectral 
type is a proxy for evolutionary status: the giant components in L&IM binaries may 
either be located on the first red giant branch (RGB), in the core He-burning phase 
(which is hardly distinguishable from the lower RGB/AGB; CH giants probably belong 
to that phase). He-shell burning early AGB (E-AGB), or on the thermally-pulsing AGB 
(TP-AGB) phase, where the s-process operates. 

Symbiotic activity is expected in the middle of this spectral sequence, because (i) at 
the left end, stars (like CH) are not luminous enough to experience a mass loss sufficient 
to power symbiotic activity; (ii) at the right end, the stars with the barium syndrome need 
not be binaries. Indeed, in TP-AGB stars, heavy-elements are synthesized in the stellar 



interior and dredged-up to the surface, so that "intrinsic Ba" (or S) stars occupy the 
rightmost boxes - 4 and 8 - of Fig. [21 and hence need not exhibit any symbiotic activity 
since they are single stars (For examples of stars belonging to box 4, see ifisll and [201). 
Such evolved giants which are nevertheless members of binary systems (like Mira Ceti 
- box 9) will of course exhibit symbiotic activity. It is noteworthy that late M giants are 
inexistent in a halo population (hence the crossed box 5), because evolutionary tracks 
are bluer as compared to higher metaiiicities. Examples of such very evolved (relatively 
warm) stars in a halo population (box 4) include CS 30322-023 Q and V Ari Q. 



3. IMPORTANT PRELIMINARY NOTIONS 

3.1. Orbital elements 

The 7 elements used to describe an orbit are the semi-major axis a, the eccentricity e, 
the orbital period P, the time of passage (at periastron for a non-circular orbit) Tq, the 
orbital inclination on the plane of the sky z, the longitude of periastron (for non-circular 
orbits) (D, and finally the position angle of the ascending node Q. The angles z, (O and Q. 
are identified in Fig. [31 

These 7 elements are also called Campbell elements, as opposed to the Thiele-Innes 
elements described below. For edge-on orbits, z = 90°. The semi-major axis may either 
refer to the relative orbit of the two components (usually denoted a, in the case of visual 
binaries), to the orbit of one component with respect to the centre of mass of the system 
(usually denoted a^, as, in the case of spectroscopic or astrometric binaries), or to the 
orbit of the photocentre with respect to the centre of mass of the system (usually denoted 
ao, in the case of astrometric binaries). For astrometric and visual binaries, the semi- 
major axis is an angular quantity (sometimes denoted a"), so that the conversion to a 
linear quantity requires the knowledge of the parallax U5: a(AU) = a" / U5. 

For computational reasons, it is often convenient to replace four of the Campbell 
elements, a,z, £0,n, by the so called Thiele-Innes elements (or constants). The apparent 
motion of a binary component in the plane of the sky (i.e., the plane locally tangent 
to the celestial sphere) is described by the cartesian coordinates (jc, y) (with x pointing 



towards the North) ulUM 



x = AX+FY 

y^BX + GY ^ ^ 



with 

X = cosE 



Y = a/ 1 — e^sin£', 

where {X, Y) are the coordinates in the true orbit and E is the eccentric anomaly related 
to the mean anomaly M by Kepler's equation 

M{t) = 27t^-^^=E-esmE. (2) 




FIGURE 3. The orbital plane and orbit of component A around the centre of mass G. z is the coordinate 
along the line of sight, pointing away from the observer. The plane of the sky is tangent to the celestial 
sphere and perpendicular to the line of sight. 



The Thiele-Innes constants A,5,F, G are related to the remaining orbital elements by 

A = a(+cos (UcosQ — sinft) sinQcos?) (3) 

B = a{+ cos CO smQ. + sin CO cos Clcosi) (4) 

F = a{— sin CO cos Q. — cos CO smQ.cosi) (5) 

G = a(— sin ft) sin fl + cos ft) cos n cos/). (6) 



3.2. Roche lobe 

The concept of the Roche lobe is central while dealing with binary stars. Let us 
consider a circular binary and a test particle corotating with this binary, such as a mass 
element on the surface of a star rotating synchronously with the orbital realiasvolution. 
The Roche lobe corresponds to the critical equipotential (in the corotating reference 
frame) surface around the two binary-system components which contains the inner 
Lagrangian point. At this point the net force acting on the corotating test mass vanishes 
(the net force is the vector addition of the two gravitational forces and of the centrifugal 
force arising in a non-inertial reference frame rotating uniformly with the binary). If the 
star swells so as to fill its Roche lobe, matter will flow through the inner Lagrangian 
point onto the companion. The radius of a sphere with the same radius as the Roche lobe 



around star A may be expressed by 



- 0.6^-2/3+in(l+^-i/3)' 

where q = Mb /Ma is the mass ratio. In the case of a star losing mass (like an AGB star), 
the extra-force responsible for the ejection of the wind must be included in the effective 
potential. It distorts the shape of the equipotentials and reduces the size of the Roche 



lobe around the mass-losing star [|26 



. To take this effect into account, a generalisation 



of Eggleton's formula is necessary Mn\- This reasoning is strictly applicable only to 
circular systems, but is often employed as an useful approximation also in the case of 
non-zero eccentricity. A more detailed discussion of the concept of Roche lobe is given 
by, e.g.. Shore [2] or Jorissen [|4i]. 



4. SPECTROSCOPY 



4.1. The method 



Spectroscopic detection of binaries relies on measuring the Doppler shifts of stellar 
spectral lines as the binary components orbit their centre of mass. Soon after the discov- 
ery of the first spectroscopic binaries in late nineteen century [28L129] this became the 
preferred detection method as it is robust and through relatively simple means it gives 
access to important physical parameters. Time-dependent spectroscopy allows direct de- 
termination of P,To,e, and CO - period is the easiest to constrain, the other parameters 
require more detailed knowledge of the shape of the radial- velocity curve. Radial ve- 
locities give no handle on Q., but convey entangled information on a and i. Joined with 
Kepler's third law, this leads to partial, but very much sought after, knowledge concern- 
ing the component masses (see Sect. 14. 3 1 and Table[5]in Sect. [8]). 

The expression for the radial- velocity of star A is 



K(A)=K(G) + dz(A)/d? = K(G) + KAiecos(o + cos{(0 + v{t))), 



(8) 



where K(G) is the radial velocity of the centre of mass of the system, z is the spatial co- 
ordinate along the line of sight (Fig. [3]), v is the true anomaly (angle between periastron 
and the true position of the star on its orbit), and [e.g. ..30.1 
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is the semi- amplitude of the radial- velocity curve of component A, with q = Mb/Ma- 
To fix the ideas, the semi-amplitude Ka associated with companions of various masses 



TABLE 1. The semi-amplitude (in 
km s^') associated with companions of various 
masses and orbital periods, for a 3 M0 primary 
star in a circular orbit with / = 90°. 
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and orbital periods is given in Table [H for a 3 primary star in a circular orbit with 
/ = 90°. 



4.2. One-dimensional observation: the uncertainty introduced by the 

inclination 

The first difficulty facing spectroscopic -binary detection comes from the fact that the 
radial velocity is a one-dimensional measurement (along the line of sight), which implies 
that the knowledge of the orbit can only be partial. In particular, what is known is the 
projected orbit on the plane of the sky, implying a degeneracy between / and ua, since 
only the combination ua sinz may be extracted from Ka- The same expression holds for 
component B if that component is visible ('SB2 systems' standing for spectroscopic 
binaries with two observable spectra). 

If the components of a binary are of approximately equal luminosities, the spectrum 
will appear as 'composite' (as further discussed in Sect. 17.31) . and a radial-velocity curve 
(Eq.[8l) will be available for both components. Since uaMa = agMs from the definition 
of the centre of mass, the mass ratio q = MB/MA maybe derived from /^^//Cg (see Eq.|9]). 
Then, Ka +Kb joined with Kepler's third law 

a^/p^ = G {Ma+Mb)/Ak^, (10) 

where a = + <5[fi is the semi-major axis of the relative orbit, gives access to {Ma + 
Mb) sin^ / (since only a sin / - not a - may be extracted from Ka +Kb). Combining the 
mass ratio and sum (multiplied by sin i) then yields Ma sin i and Mb sin i for the case 
of SB2 binaries. A summary of what may be known about the masses for spectroscopic 
binaries is provided by Table |5] in Sect. [8l 

If there is an independent way to derive the orbital inclination i (the system is visual, 
astrometric or eclipsing; in the latter case i is close to 90°), then (and only then) 
may the individual masses be derived (see Table [5]). A textbook case [among many 
others; see the review in 'sT] is provided by the S star HD 35155 : this star is a 
spectroscopic binary with only one spectrum visible in the optical region [32J, but an 
International Ultraviolet Explorer (lUE) spectrum reveals ultraviolet emission lines tied 
to the companion which thus gives access to the mass ratio. Eclipses have been 



observed in lUE spectra and optical photometry [|34L 1351]. which implies an inclination 



close to 90°. With all these data at hand, the masses inferred for the giant and its 
companion are 1.3 - 1.8 M© and 0.45 - 0.6 M©, respectively [34]. 

If one is only interested in the distribution of the masses (rather than in the masses 
of individual objects) for a sample of SB2 binaries with their orbital planes inclined 
randomly on the sky, then statistical techniques may be used, as further discussed under 
the next item. 



4.3. Only one spectrum is observable: the mass function 

It must be stressed that the semi-major axis a = UA+ag of the relative orbit cannot be 
derived if only one spectrum is observable (for instance in the case of a faint companion, 
on the lower main sequence or a WD). Hence Kepler's third law, which involves a, is 
not applicable. Instead, for these binaries with only one observable spectrum (say A), 
only a quantity called the mass function and denoted f{M), can be derived: 

, , (Mr sin kIp , 

= 1.036 X 10"^ if](kms-i)i'(d) (1-^2)3/2^0 (12) 

Still, there is a way to get better constraints on the masses if specific conditions 
are met, namely Mb « Ma (exoplanet companion) with Ma known from the mass- 
luminosity relationship for main-sequence stars (if applicable). In those circumstances, 
Jorissen et al. [36jl have shown that, for a sample of stars, it is possible to extract 
the distribution of Mg from the observed distribution 4>(y) of Y{Mb) = Mssini. This 
follows from 

A/n sin"^ / 

f{M) ^ ^ , (since Mb «Ma). 

A 

As Ma is known, Y{Mb) may be derived from 7 (Mg) = f{M) 

In the more general case of a stellar rather than planetary companion, and with Ma 
known (as above), the distribution of 2 = ^ (1 +^)^2/-^, where q = Mb/Ma, may be 
obtained rather than that of Mb [40] . Indeed, we have f{M) =MAq^ (l+q)^^ sin^ i = 

Ma sin^ /. Hence, Y{ff) = Qsini = f{M) ^/^M^ is available from the observations. 
Then, the sought distributions ^(Mg) or ^{Q) obey the relation: 

noo /'OO 

4>(y)=/ ^(Mfi)n(y|Mg)dMfi or 4>(y)=/ ^(2) n(y|e) dg. 

Jo Jo 

The kernel II{Y\Mb) corresponds to the conditional probability of observing the value 
y given Mb- Under the assumption that the orbits are oriented at random in space, the 



^ Alternative methods have been proposed by Mazeh and Goldberg IstIi . Zucker and Mazeh jSSll and 
Tabachnik and Tremaine ll39ll 



inclination angle i distributes as sinz : n(z) di = sinz dz and the expression for the kernel 
is obtained from: 

n(Mesinz')A/gcosz dz = sinz dz. 
Eliminating the inclination in the above relation yields 

and 

*(!-) = ^/^(Mri-j^-^j^jl^dM.. (14) 

This integral equation must be solved for ^(Mg). It can be reduced to Abel's integral 
equation by the substitutions 114 ill 

= 1/77 and M| = l/t. (15) 
With these substitutions, Eq.([T4l) becomes 

where 

HV)^^{^) and VA(0 ^ (17) 
The formal solution of Abel's equation (Eq.[T6l) is given by 

n Jo dri {t-riy/^ n y/t 

where ^{0) = limy^.c4>(y) = 0. 

It is difficult to implement numerically, since it requires the differentiation of the 
observed frequency distribution ^{Y). Unless the observations are of high precision, 
it is well known that this process can lead to misleading results. Two approaches are 
possible to overcome that difficulty: 

• The observed frequency distribution is smoothed in an optimal way before being 
used in Eq.[T8] 14211 . The solution ^{t) is then computed numerically using standard 
differentiation and integration schemes. 

• The Lucy-Richardson inversion algorithm is applied to Abel's integral equation 
110,43,44]. 



Both methods have been used by Jorissen et al. Bql to derive the distribution of exo- 
planet masses. 
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FIGURE 4. AAVSO light curve (filled dots) for the Mira variable R CMi during three successive cycles, 
compared to its radial-velocity curve (solid line and open dots, from ll32ll ). Radial-velocity measurements 
were not necessarily obtained at the indicated dates, as they were all folded onto the radial-velocity 
solution. Note that the light- and velocity-curves are phase-shifted: maximum light occurs when the 
velocity curve crosses the centre-of-mass velocity (represented by the horizontal dashed line). (From Q]) 

4.4. Intrinsic velocity variations 

4.4.1. Pseudo-orbital variations and Roche lobe radius 

The pulsation of the atmosphere of Mira, Cepheid or RV Tau variables causes intrinsic 
velocity variations and makes binary detection using radial velocities difficult or even 
impossible. Often the radial- velocity variations associated with the pulsations are with 
the same period as the photometric variations (although they are phase-shifted), so that 
the intrinsic nature of these radial- velocity variations is very clear. A very spectacular 
example is provided by the comparison of the light curve of the 337 d Mira variable 
R CMi with its radial- velocity curve, as shown on Fig. IH The radial- velocity semi- 
amplitude is as large as 8 km s^^. Clearly, in such a case, a companion could only be 
detected if it would yield orbital variations of at least the same order (compare with the 
data of Table [H to see which kind of companions would be detectable). 

Hinkle et al. have shown that the radial-velocity semi- amplitude associated with 
Mira and semi-regular pulsators correlates well with the visual amplitude and with the 
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FIGURE 5. The radial-velocity curve of the RV Tau variable IRAS 08544-4431: the valuations with a 
large amplitude are caused by the orbital motion, whereas the small-amplitude variations are caused by 
the RV Tau-like pulsations with period 90 d, as shown on the bottom panel, displaying the residuals with 
respect to the orbital solution. (From Ii46l1 ) 



pulsation period, reaching 15 km in the most extreme cases. 

The binary post-AGB star IRAS 08544-4431, hosting a RV Tau variable (a class of 
luminous variables in the Cepheid instability strip, characterized by alternating deep and 
shallow minima), is a rare case where intrinsic and orbital radial-velocity variations are 
superimposed (Maas et al. 2003 i46|l present in fact quite a number of similar cases). 
Here the orbital period is 499 d and the orbital semi-amplitude is 8 km s ^ to be 
compared to 4 km s^^ for the pulsations of period 90 d (Fig. [5]). Hence the two kinds of 
variations may be separated, but at the expense of a dense observational coverage. 

Sometimes, however, the situation is not as clear as for the two cases discussed above. 
The 225.9 d Mira variable S UMa exhibits radial velocity variations that could at first 
be interpreted as orbital motion with an apparent period of 576 d and an eccentricity of 
0.29 1I32I1 . based on a rather scarcely sampled data set (Fig. (6]). However, the last two data 
points deviate markedly from the solution based on earlier data points, casting doubts 
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FIGURE 6. Left panel: The radial-velocity curve of the 222 d-period Mira variable S UMa. A solution 
with a period P — 576 d and an eccentricity of 0.29 provides a good match to all the data points but the last 
two. Since the 1/576 d^' frequency is close to the 1/222.0 - 1/365.25 alias of the 1 -yr frequency and of 
the pulsation frequency, the intrinsic origin of the radial-velocity variations becomes the best explanation. 
From ll47ll . Right panel: The radial-velocity data of S UMa phased with the 222 d pulsation period. The 
solid line is an (unsatisfactory) keplerian solution of eccentricity 0.29. 



on that solution. A period analysis of this data set 14711 reveals that the 576 d period is 
probably an alias of the pulsation period (present in the radial-velocity data as 222.0 d) 
and of 1 yr (1/222.0 - 1/365.25 = 1/566.0), thus strongly reinforcing the suspicion that 
the radial-velocity variations have an intrinsic origin. 

la any case, periods as short as a few hundred days can in no way be associated with 
an orbital motion involving a star as large as a Mira, since the Mira would then fill its 
Roche lobe, and the system would exhibit strong signatures of mass transfer (like for 
instance an X-ray flux, or strong emission lines). To hold within a binary system, the 
radius i? of a Mira variable derived from the relationships [|48ll 

logi?Mira = (logP 0.9 logM-F 2.07)/ 1.9 fundamental mode, 

(19) 

logi?Mira = (logP -f 0.5 logM+ 1.40)/ 1.5 first overtone mode, 
where P is expressed in days, and R and M in solar units, must be smaller than the critical 



Roche lobe radius Rr expressed by Eq.r7][|25l.l27[]. The orbital period for which a Mira of 



a given pulsation period (in either the fundamental or first overtone mode) fills its Roche 
lobe is displayed in Figure |7j It is derived from the above formulae ([7]) and (fT9l ), and 
from Kepler's third law (assuming typical masses for the Mira star and its companion, 
as indicated on the figure). 

The orbital periods allowed by this criterion are quite long (more than 1000 d), and are 
always much longer than the pulsation periods. The corresponding radial- velocity semi- 
amplitudes K (Fig. |7]) may be computed from Eq. [9l This quantity must be compared 
with the intrinsic radial-velocity variations of Mira variables (due to pulsations) which 
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FIGURE 7. Bottom panel: Orbital periods versus pulsation periods for Miras in semi-detached sys- 
tems, derived from the period-radius relations (Eq. [19] ) for the fundamental and first overtone modes. 
Detached systems are located above the corresponding lines. The filled circles labeled MO and M7 indicate 
the critical periods corresponding to the median radii of non-Mira MO and M7 giants |49], respectively 
(adopting masses of 1.0 and 0.75 Mq for the giant and its companion, respectively). Top panel: Orbital 
radial-velocity semi-amplitude versus pulsation periods (thick lines: fundamental mode; thin lines: first 
overtone) for Miras in circular semi-detached systems with sin; = 1 (where ; is the orbital inclination). 
Binary systems with larger semi-amplitudes would involve Roche-lobe-filling AGB stars. (From |j4|]) 



can be as high as 15 km s ^ ll45ll . The detection of spectroscopic binaries among Mira 
variables thus appears to be almost impossible. 



4.4.2. Wood's sequence D and long secondary periods 

Variable stars with 'long secondary periods' (LSP, characterized by periods between 
200 and 4000 days, i.e., a factor 5 to 15 longer than the primary period) and with V- 
band amplitudes up to 1 mag have been known for decades [50, 51], but the interest for 
it has been renewed since Wood [1521 . 15311 showed that these secondary periods follow a 
period-luminosity (PL) relation (the so-called 'sequence D') in the PL diagram of long- 
period variables (LPVs) in the LMC. Besides the expected sequence corresponding to 
the fundamental pulsation mode of Miras (the so-called 'sequence C) and its higher 
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FIGURE 8. Period-luminosity relations for red giant pulsators discovered by the MACHO project in 
the Large Magellanic Cloud. Ellipsoidal and contact eclipsing binaries are shown as plus signs. The solid 
lines show the critical orbital periods for systems of unit mass ratio filling their Roche lobe, with primary 
masses of 0.85 and 2.5 Mq. (From Ii54l ) 



harmonics (sequences A, A' and B, mainly populated by semi-regular variables), a 
sequence very clearly appeared at longer periods (the so-called 'sequence D'; Fig. [8]), 
involving about 25% of the LPVs in the LMC. Various hypotheses have been proposed to 
explain the origin of the LSP variability: rotation of a spotted star, episodic dust ejection 
or obscuration, radial or non-radial pulsations, stellar or substellar companions. Wood 
et al. 1521 ] suggested that stars on sequence D are components of semi-detached binary 
systems. But is it possible that 25% of all the LPVs in the LMC are member of binary 
systems, let alone semi-detached systems? 

Derekas et al. and Soszyhski et al. fs^] noticed that the so-called 'sequence E', 
containing eclipsing and ellipsoidal binaries, is the extension of sequence D for fainter 
(i.e. non LPV) stars. Soszyhski [|56ll obtained another clear result calling for the role 
of binarity in the LSP phenomenon: by considering the quantity \3Rm — Bm (where 
Rm and Bm are the MACHO magnitudes in the red and blue bands, respectively), the 
LSP phenomenon should disappear since Derekas et al. [54] found that LSP stars have a 
Bm/Rm amplitude ratio of 1.3 on average, as compared to 1.1 for ellipsoidal variability, 
which should be grey. Soszyhski [|56ll found that in 5% of the stars belonging to se- 
quence D, a double-humped structure with a period equal to the LSP one is observable 




FIGURE 9. The structure of the flow (represented by its velocity field) from a mass-losing AGB star 
(located at x = 0.33, y = 0) around the accreting star (located atx = —0.66,y = 0) in a frame rotating with 
the system. The double-spiral structure is clearly visible in this simulation of an adiabatic flow. (From 
Hi) 



in the 1.3Rm — Bm residual lightcurve, but with a slight phase lag with respect to the 
LSP lightcurve. This finding gives credit to the binary hypothesis, and the phase shift 
even hints at dust obscuration as being responsible for the double-humped lightcurve, 
since hydrodynamical simulations of the mass transfer involving a mass-losing AGB 
star reveals that the accretion column forms a spiral arm behind the companion [Fig. |9j 



571 |58L |59, 60, 61]. Therefore, the eclipse of the star by this matter trailing behind the 



companion occurs slightly before the geometric conjunction. 

This will be discussed further below in Sect. 15. 3[ where we will show that such dust 
obscuration events are rather common among binaries involving mass-losing giants. The 
dust cloud present in the LSP stars is not very conspicuous, though, since mid-infrared 
colors of LSPs and non-LSPs are similar [62, [63l] and there are no LSPs showing the 
large [60]-[25] fim IRAS excesses exhibited by some R Coronae Borealis stars, a class 
of stars known for forming dust in large quantities. 

Another piece of evidence in favour of binarity was presented by Derekas et al. [|54ll . 
Soszyhski et al. llssll . Soszyhski [f56ll and Jorissen et al. i64ll who showed that Wood's 
sequence D closely matches the sequence expected for binary stars involving giant stars 
filling their Roche lobe. This is clearly illustrated by the solid lines in Fig. [8] and by 
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FIGURE 10. The orbital-period as a function of (absolute) K magnitude for galactic M giants in binary 
systems. The ordinate axis corresponds to the K magnitude that the M giants would have if they were put 
at the distance of the LMC (see text). The dashed line labelled 'D' corresponds to Wood's sequence of 
long secondary periods. The solid and dotted lines correspond to the limit on the orbital period imposed 
by the Roche lobe radius (for M\ = 1.3 M0 and M2 = 0.6 M©), for stellar radii comprised between 40 
and 120 R©. (From [64])- 



Fig. [TOl which is a diagram orbital-period - luminosity for the exhaustive sample of 
spectroscopic binary M giants from Jorissen et al. [64]. The ordinate axis of Fig. [10] 
corresponds to the K magnitude that these M giants would have if they were at the 
distance of the LMC, i.e., K{IMC) = Mk + 18.50, where the McNamara et al. 
LMC distance modulus has been adopted, and the absolute K magnitude is derived from 
the 2MASS value combined with the distance from the maximum-likelihood estimator 
of Famaey et al. [[66|] . It is very clear that Wood's sequence D does indeed come close 
to the upper envelope of the region occupied by the galactic binary M giants, which is 
defined by the condition that they hold within their Roche lobe. 

If Wood's sequence D is indeed related to LPVs filling their Roche lobe, the remaining 
question is: how come that there are so many? Soszyhski [|56ll suggests that variability 
along sequence D originates in giants with substellar companions, the latter being 
supposedly much more numerous than the 15% of spectroscopic binaries (not even 
restricting to semi-detached systems) observed among M giants [|67|] . Although the mass 
functions of the orbital solutions obtained by Wood et al. [|63n and Hinkle et al. ll68n 
are indeed compatible with substellar companions, the eccentricities and longitudes of 
periastron derived for 5 among the 6 orbits presented by [68] are surprisingly similar 
and may cast doubts on the orbital origin of the observed variations. 



5. PHOTOMETRY 



As stated in the previous section on spectroscopic binaries, not many spectroscopic bina- 
ries are known among Mira variables, because the orbital radial-velocity semi- amplitude 
^orb corresponding to a long-period binary is generally smaller than ^puis- Other meth- 
ods must thus be used to detect binaries involving Mira variables! Fortunately, there are 
specific methods using photometry to do so. 



5.1. Miras with flat-bottom light curves 

Miras exhibiting a light curve with a flat bottom and a small-amplitude, despite their 
Mira-like period, are good candidates for being binaries with a faint companion, which 
dominates the system light around minimum light and is responsible for the flat bottom. 
A good example thereof is Z Tau {P ~ 450 d), whose AAVSO (see Sect. (9]) light curve 
exhibits a flat bottom at V = 14. 



5.2. Ellipsoidal or eclipsing variables 

Ellipsoidal variables are characterized by a light cycle which is exactly half the orbital 
period, caused by the tidal deformation of the star nearly filling its Roche lobe. The 
detection of ellipsoidal variations or of eclipses has often been the first evidence of 
binarity for many systems. It offers an easy way to infer the binary period [see 69, for 
the case of symbiotic stars]. Ellipsoidal variations should have the same amplitude in all 
photometric bands, as it is a purely geometrical effect ff^ . 



5.3. Dust obscuration and circumbinary disks 

Hydrodynamical simulations of a wind-driven accretion flow in binary systems |5 
59I 60, 61, 71] predict the formation of a spiral pattern corresponding to the accretion 



wake bended by the Coriolis force (Fig. [9l). This prediction has been nicely confirmed 
by the direct imaging in visible light of the proto-planetary nebula AFGL 3068, using 
the ACS camera on board the HST, which reveals a spiral pattern winding several times 
around the central star [72]. Periodic obscuration of the central star by dust clumps 
trapped in the accretion structure - be it spiral wave, disk around the companion, or 
a circumbinary disk - seems common in binary systems involving a mass-losing giant 
star. This phenomenon very likely plays a role in the photometric variability observed in 
the following classes of stars: 
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FIGURE 11. The V band light curve of V Hya between October 1961 and July 1996 from the AAVSO 
archives. These data show a more or less regular variability with a period of 533 d and a peak-to-peak 
variation of 2 magnitudes at V, plus a longer-term variation with a period of about 6500 days (17-18 
years) with deep minima (5 to 6 magnitudes). (From 1,75.1 ) 



Binary systems with evidence for dust obscuration 



Stellar class 


Examples 


References 


LSP Wood's sequence D 


many 




short Porh Ba or S stars 


HD 35155, HD 46407, HD 121447 


[34,73,74] 


AGB C star 


VHya 


[75J 


Post-AGB 


HR 4049 


\M 


Supergiants 


e Aur 


[77] 



The case is especially clear for the C star V Hya (Fig.[lTI). Quoting Knapp et al. 117 5h : 

The morphology of the 17 y variation resembles that of an eclipsing binary, but with an 
eclipse duration which is far longer than can be produced by a stellar companion and 
of an amplitude which shows that essentially the entire stellar photosphere is occulted. 
We suggest that the regular long-period dimming of V Hya is due to a thick dust cloud 
orbiting the star and attached to a binary companion. The orbital period is inferred from 
the light curve displayed on Fig. [HI and there is no spectroscopic confirmation so far 





disk erosion due to 
radiation pressure 



FIGURE 12. The V light curve of the post-AGB star HR 4049 (upper left panel). The modulation is 
caused by the line of sight going through the dusty circumbinary disk once per orbit, according to the 
model depicted in the bottom panel. The upper right panel compares the amplitude of the variations in 
various photometric bands with the extinction law in the interstellar medium. (From |,76,] and I.8Q1I ) 



(difficult because of the long period, hence small semi-amplitude!). 

The 27-yr period observed in the FO supergiant e Aur [77] may be a similar case, 
with the long duration of the eclipse showing that the secondary cannot be a star but 
must be a large (510 AU) dark body, which observations strongly suggest is a dust disk 
orbiting a companion [IttI , ItsIi . Eclipses best explained by dust clouds attached to a 
binary companion are also seen in a small number of planetary nebula nuclei, e.g. in 
NGC 2346 [79]. Finally, the prototypical post-AGB star HR 4049 exhibits eclipses by a 
dust disk, which has been imaged directly by VLTI [Fig.[l2ll76|,|80|,|8i|]. The geometry 
of the binary system (bottom panel of Fig. [T2l) is such that the line of sight goes through 
the circumbinary disk once per orbital cycle, thus the resulting light curve has only 
one minimum per cycle (upper left panel of Fig. [T2l) rather than two as in ellipsoidal 



variables. The color dependence of the amplitude of the variations is the same as the 
interstellar extinction law (upper right panel of Fig. [T2l), thus confirming that the optical 
variability is caused by scattering on dust grains of similar size as in the interstellar 
medium. 

Such circumbinary disks are almost a defining property of binary post-AGB stars, 
since the presence of a 21 /im excess, caused by the presence of cold dust in a circumbi- 
nary disk extending far away from the heat source (the post-AGB star), may be used for 
identifying binary post-AGB stars [fl4l [82i.[83ll. 



The M4III -I- A system SS Lep (= HD 41511) shares with post-AGB stars the presence 
of a circumbinary disk 1841 ISSD, fed by non conservative RLOF. Using interferometry, 
Verhoelst et al. [85] have indeed shown that the M giant fills its Roche lobe, and that 
the A type companion has a radius about 10 times larger than that expected for a main 
sequence star. It has very likely swollen as a result of accretion. The system HR 1129 
(G2Ib -I- B7) shows as well an infrared signature of dust most likely associated with a 
circumbinary disk H^. 

The lightcurves of the Ba star HD 46407 and the S star HD 35155 (see Fig. [H) are 
a bit more puzzling. Barium stars are post-mass-transfer objects with WD companions 
[ItI. [stL 88 1 which are too faint to yield detectable eclipses in the visual. And yet, the 
reality of the eclipses in the short-period Ba star HD 46407 (P = 457 d) cannot be 
doubted, since the eclipses have been observed during different cycles, though with 
variable depths and phase lags [|73L l741. A very important clue as to the origin of these 
eclipses comes from the fact that the deepest eclipse (observed around JD 2447 400) has 
occurred just before an episode of secular obscuration (left panel of Fig. [141) . The fact 
that the color dependency of the amplitude for both the eclipse and secular variations 
follows a law (right panel of Fig. [141). reminiscent of Rayleigh scattering, at least in 
the Stromgren u and v bands, is a further indication that the phenomenon is caused by 
(small) dust grains. In b and y, it is closer to the ISM color absorption law. Despite the 
fact that barium stars host K giants which do not lose large amounts of mass, the secular 
obscuration event calls for the presence of recent dust in the system, rather than for the 
remnant of the spiral arm dating back to the epoch where the current WD companion 
was a mass-losing AGB star. 

A similar eclipsing behaviour has been reported for the S star HD 35155 with 640 d 
period by 04I1 . Despite the fact that Adelman 18911 did not confirm the occurrence of the 
eclipse, the two lUE spectra of HD 35155 obtained by chance at the time of eclipse have 
the smallest flux among the 10 available lUE _spectra, though neither the continuum nor 
the emission lines disappear completely 11331 13511. 

HD 121447 is the coolest barium star (K7III) with the second shortest orbital period 
(186 d). Ellipsoidal variability was first suspected for this star [90], but the absence 
of synchronous rotation and a strong color dependence of the light curve [89] do not 
support this conclusion, so that light variations caused by scattering on dust formed by 
the compression of the red-giant wind in the accretion process, becomes an attractive 
possibility. 

To summarize, disks (often circumbinary) are observed in many different classes of 
binary stars involving a mass-losing star, and seem to be very common in such circum- 
stances. They play an important role in the evolution of such systems: 
- they control the evolution of the eccentricity through tidal effects and angular momen- 
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FIGURE 13. The light curves phased with the orbital period for the barium star HD 46407 (left panel; 
P = 457 d) and the S star HD 35155 (right panel; P = 640 d) which are among the systems with the 
shortest periods in these classes. Data for HD 46407 [|73iF74ll and HD 35155 [34] in the Stromgren system 
from the Long-Term Photometry of Variables program ll9lll . The small arrows in the right panel identify 
the phases corresponding to the eclipse (0.17) and transit (0.65) of the companion. HD 35155 exhibits as 
well intrinsic photometric variations. 
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FIGURE 14. Left panel: The long-term light curve (in the Stromgrem bands uvby, from bottom to top) 
of HD 46407 with respect to the average of the comparison stars (left), and the difference between the two 
comparison stars (right). The solid line roughly illustrates the secular change in the average brightness. 
See Jorissen et al. 07411 for details. Right panel: The amplitude of the secular and eclipse variations in the 
M, V and b bands normalized by the amplitude in the y band, and comparison with the interstellar-medium 
extinction law and the Rayleigh law. 



turn exchange with the orbit ll92ll ; 

- they trigger a dust/gas segregation, so that the abundance pattern observed at the surface 
of the post-AGB star is shaped by the condensation temperature of the various elements 
[|93|] . It is believed that this physical segregation operates through re-accretion by the 
post-AGB star of gas depleted in refractory elements which stayed in the dust phase. 
It leads to very Fe-depleted post-AGB atmospheres (down to [Fe/H] = —4.8) [94]. In- 
triguingly, a similar scenario has recently been proposed n95\\ as the possible cause of 
the ultra low metallicities observed in several stars from the Hamburg/ESO survey [|96|l . 
like the record holder HE 1327-2326 with [Fe/H] = —5.6 11951 l97ll. This scenario would 
require these stars to be binaries. There is no definite evidence thereof so far, but the 
available radial- velocity data is not necessarily conclusive [j95l]. 



6. ASTROMETRY 
6.1. The basics 

Very long-period systems are hard to detect with spectroscopy or photometry. This is 
where astrometry comes to rescue. Astrometric binaries differ from the more traditional 
visual binaries in the sense that for visual binaries, both components are observed so that 
it is the relative orbital motion (of one component with respect to the other) which is 
readily detected. For astrometric binaries instead, it is the absolute (non-linear) motion 
of one component on the sky which is detected. Space astrometry, which started with 
the Hipparcos satellite [98] launched in 1989, has reached the milliarcsecond (mas) 
accuracy level for stars down to about V = 10, and the coming Gaia satellite (to be 
launched in 201 1) should improve this accuracy by a factor of 100. Many achievements 
have already been made in the field of binary stars by Hipparcos, as reviewed by 
Ferryman |99|l . and many more must be expected from Gaia. We will discuss just a 
few here, illustrating the potential of astrometry in the discovery of specific kinds of 
binaries [A very extensive review of this potential is presented in 99]. For instance, 
since astrometric methods are independent of spectral types (as opposed to spectroscopic 
methods, since the detection of spectroscopic binaries is dependent upon the possibillity 
to follow accurately the variations of the position of spectral lines), astrometry offers 
a way to derive unbiased frequencies of binaries of different spectral types [Table [31 
QOO]. Su bstellar compa nions may also be detected thanks to astrometry, if accurate 
enough lIlOll Il02l [l03il . Astrometry offers as well good prospects to detect binaries 
involving Mira variables (which are especially difficult to find with spectroscopy), 
provided however that spots or an inhomogeneous surface brightness or asymmetries 
in the shape of the stellar disk of these very extended stars do not cause var iations of the 



position of the photocentre that confuse the parallactic and orbital motion [11041 . 1 1051] . 

The principle of the astrometric measurement by Hipparcos and Gaia needs to be ex- 
plained first. In order to achieve a good accuracy overall on the celestial sphere, Hippar- 
cos had two fields of view, 0.8 degree square each, widely separated on the sky (by about 
58 degrees of arc) and superimposed on the detector. In the case of Gaia, the numbers 
will be 0.37 degree square and 106.5 degrees, respectively (left panel of Fig. [T5l) . Since 
the satellite is spinning, this signal from stars crossing the instrument fields of view is 




FIGURE 15. Left panel: The scanning law of Gaia. (From the ESA-Gaia web site) Right panel: 

The modulating grid in the Hipparcos focal plane (top), for two different scanning directions A and B, 
producing the modulated signals displayed on the bottom panel. The phase p for the source S is measured 
differentially with respect to that of the reference point R, modulo the grid step. In reality, the grid extends 
over the images of S and R in the focal plane, although this was not displayed for clarity. (From 1110^ ) 
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FIGURE 16. Definition of the abscissa offset v for a predicted position {a,,,, 8m), with respect to the 
observed position somewhere along the thick line perpendicular t o the reference great circle labeled 
'RGC. Remember that the measurement is one-dimensional. (From llOTll ) 



modulated by a one-dimensional grid in the focal surface (right panel of Fig. [15]). As 
the satellite scanned the sky in a complex series of precessing great circles, maintaining 
a constant inclination to the Sun's direction, a continuous pattern of one-dimensional 
measurements was built up. In order words, it is the 'abscissa' along the great circle 
which is being measured, the position of the star along the perpendicular to the ref- 
erence great circle being unknown. These one-dimensional abscissa measurements are 
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FIGURE 17. Left panel: A solution fitting the asti-ometric motion of HIP 85749 with a quasi-parabolic 
orbit (e = 0.998) aligned with the line of sight. Right panel: The orbital solution with e = 0.2, consistent 
with the spectroscopic orbital elements. Both solutions have similar goodness-of-fit values. The small 
line segments connect the predicted position with the observed great-circle abscissa (see Fig. [15). (From 

Coi) 



then confronted to the expected positions of the star on the sky as a function of time, 
for a given reference position, parallax and proper motion for that star. The difference 
between these one-dimensional positions is called 'abscissa residual', and is noted v on 
Fig. [m [107]. These abscissa residuals are then combined with the variance-covariance 
matrix of the obse rvational errors to yield a 'x} value characterizing the quality of the 



solution [[lOTLllQSl] 



Right away, two difficulties become apparent in the context of binary stars: 
(i) It is the position of the photocentre of the system which is measured, reducing to 
the position of the brightest component if the difference in brightness is large. It may 



be shown [|23|] that the relation between the semi-major axis qq of the orbit of the 



photocentre around the centre-of-mass and semi-major axis a of the relative orbit writes 

aQ = a{K-^), (20) 

where K = A/2/ (Mi + M2) and /3 = 1 / ( 1 + loP-'^^^^^ where Am = m2 - mi , so that when 
component 2 is much fainter than component 1 (i.e.. Am — > 00), /3 = and aQ = aKh then 
the semi-major axis of the orbit of component 1 around the centre-of-mass. Therefore 
the access to the masses is not straightforward from astrometric orbits. A summary of 
what may be known about the masses for the various kinds of binaries (and combinations 
thereof) is listed in Table [5] in Sect. [8] 

(ii) The one-dimensional nature of the astrometric (Hipparcos or Gaia) measurement 
imposes further limitations on the derivability of the orbital elements. As discussed 
by Pourbaix [il09ll . two-dimensional observations make it possible to draw an orbit 
projected on the plane of the sky, and to derive the areal constant Y' of Kepler's second 
law corresponding to the motion of the star on that projected orbit. Since the law of 
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FIGURE 18. Solid line: The distribution of orbital inclinations (from the Thiele-Innes set o f orbital 
elements) for the 1385 orbits from the Ninth Catalogue of Spectroscopic Binary Orbits [|l 1 ill with an 
Hipparcos entry. Many of those astrometric orbital solutions are spurious, with inclinations c lose to 90°. 
Dashed hne; The same for the solutions passing all consistency tests defined by Jancart et al. ill lOll . 



areas holds in both the true and projected orbits, one may write T' = Fcos/, where F 
is the areal constant in the true orbit, which is related to the orbital elements through 
the relation F = a^(l — e^Y^'^lTz/P and may thus be derived from these. The ratio F'/F 
then yields the inclination. With one-dimensional measurements, F', and therefore /, 
cannot be derived accurately. In those circumstances, the derived F' value is often close 
to 0, so that i close to 90° follows (edge-on orbits). To accomodate the limited arc 
span on the sky with this spurious edge-on orbit, there is no other possibility than an 
eccentricity very close to unity (quasi-parabolic orbit), a very large semi-major axis and 
a longitude of periastron close to 90° (apsidal line aligned with the line of sight). This is 
illustrated in Fig.[l7l which compares two possible solutions (having similar goodness- 
of-fit values) fitting the astrometric motion of HIP 85749. The left panel is a solution 
with a quasi-parabolic orbit (e = 0.998), an inclination close to 90° and a large semi- 
major axis, whereas the right panel displays the solution with e = 0.2 consistent with the 
spectroscopic orbital elements. 

Even when spectroscopic orbital elements are known, spurious astr omet ric orbits with 
an inclination close to 90° may still arise, as found by Jancart et al. jlloL see Fig. [TSl. 
These authors have designed statistical tests to identify those spurious solutions. A more 
detailed description of these astrometric methods to detect binaries is given in the next 
section. 



6.2. Detecting binaries from astrometric data 



A tailored reprocessing of the H ipparcos (and in the future, Gaia) Intermediate Astro- 
metric Data [hereafter IAD; 110711 makes it possible t o look for or b ital s ignatures in the 
astrometric motion, following the method outlined in jlOSll . [1112D, 111311 and applied to 



barium stars in and [ 114 1 



The basic idea is to quantify the likelihood of the fit of the Hippar cos I AD with an 
orbital model. For that purpose, Pourbaix and Arenou [115] [see also llio! ] introduced 



several statistical indicators to decide whether to keep or to dis card an orbital solution. 



The relevant criteria are as follows (we keep the notation from lIllOll ): 



The addition of 4 supplementary parameters (the four Thiele-Innes orbital con- 
stants; see Eq. \T\ and Sect. 13.11) describing the orbital motion should result in a 
statistically significant decrease of the for the fit of the A^^ IAD with an orbital 
model with 9 free parameters (Xj-^^ compared to a fit with a single-star solution 
with 5 free parameters (xh ^ parameters being the the positions a, 5, the proper 
motions /i^, jig and the parallax 07). This criterion is expressed by an F-test: 

Pr2 = Pr[F{4,N-9)>F], (21) 

where 

F = ^ (22) 
4 Xt 

Pr2 is thus the first kind risk associated with the rejection of the null hypothesis: 
''there is no orbital wobble present in the data". 

Getting a substantial reduction of the x^ with the Thiele-Innes model does not nec- 
essarily imply that the four Thiele-Innes constants A,B,F, G are significantly dif- 
ferent from 0. The first kind risk associated with the rejection of the null hypothesis 
"the orbital semi-major axis is equal to zero" may be expressed as 

Pn=Pr[x\^)>XABFG]^ (23) 

where 

Xabfg = Xs^Xt (24) 

and X^(4) is a random variable following a X^ distribution with four degrees of 
freedom. 

For the orbital solution to be a significant one, its parameters should not be strongly 
correlated with the other astrometric parameters (e.g., the proper motion). In other 
words, the covariance matrix V of the astrometric solution should be dominated by 
its di agonal terms, as measured by the efficiency £ of the matrix being close to 1 
lllldl . The efficiency is expressed by 



k^i k (25) 



nf_,v. 



k=rkk 

where and are respectively the eigenvalues and the diagonal terms of the 
covariance matrix V. 

In other words, an orbital solution will be most significant if Pr2 and Pr^, are small 
and £ large. With the above notations, the requirement for a star to qualify as a binary 
was defined as 

a = {Pr2+Pr3)/e < 0.02, (26) 




FIGURE 19. Left panel: Barium stars (previously known to be SBs) flagged as astrometric binaries 
by the algorithm are represented by black symbols. Arrows denote stars with only a lower limit available 
on the period. Right panel: Fraction of stars flagged as binaries (shaded histogram) compared to total 
numbe r of s tars, as a function of parallax. As expected, the detection rate becomes high for (D > 5 mas. 
(From El) 



where the threshold value of 0.02 has been chosen to minimize false detect ions, as 
derived from the application of the method to barium stars [see below and 111411 . 

Hipparcos data are, however, seldom precise enough to derive the orbital elements 
from scratch. Therefore, when a spectroscopic orbit is available beforehand, it is ad- 
vantageous to import e, P, Tq fror n the spec t rosc opic orbit and to derive the remaining 



astrometric ele men ts [as done in llOSLIllOl . 111211 . The above scheme has been used by 



Jancart et al. ilTToli to search for an orbital signature in the astrometric motion of bi- 
nary systems belonging to the Ninth Catalogue of Orbits of Spectroscopic Binaries [S^g; 
inn , and gives best results for systems with orbital periods in the range 100 - 3000 d 
and parallaxes in excess of 5 mas. 

If a spectroscopic orbit is not available, trial {e,P,To) triplets scanning a regular 
grid (with 10 < P{d) < 5000 imposed by the Hipparcos scanning law and the mission 
duration) may be used. The quality factor a (Eq. [26l) is then computed for each trial 
(e, P, Tq) triplet, and if there exist triplets yielding a < 0.02, the star is flagged as a binary 
(In fact, scanning only P with e set to zero allows to save considerable computing time 
and would only miss very ecc entri c binaries). To test its success rate, this method has 
been applied by Jorissen et al. on a sample of barium stars. Barium stars constitute 
an ideal s ample to test this algorithm, because they are all members of binary systems 
[Ill7llll8[|. wi th periods ranging from about 100 d to more than 6000 d. The catalogue of 
Lii et al. llll9ll contains 163 bona fide barium stars with an Hipparcos entry (excluding 
the supergiants HD 65699 and HD 206778 = £ Peg). When C7 > 5 mas and 100 < 
P(d) < 4000, the (astrometric) binary detection rate is close to 100%, i.e., the astrometric 
method recovers all known spectroscopic binaries (Fig.[T9l). When considering the whole 
sample, the detection rate falls to 22% (= 36/163) because many barium stars have small 
parallaxes or very long periods. Astrometric orbits with P > 4000 d (> 1 1 yrs) can 
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FIGURE 20. The a statistics (see Eq. l26l l as a function of the trial orbital period (assuming e = 0) for 
the P = 629 d spectroscopic binary HIP 30338 (left panel). The vertical arrow indicates the spectroscopic 
period, which falls in the region of small a. By comparison, a period of ^ 500 d may be inferred for 
HIP 90316 (right panel). The vertical dashed lines represent multiple, or integer fractions, of 1 yr At 
those periods, there is a strong correlation between the parallactic and orbital signals, which degrades the 
a statistics and makes binaries difficult to find. 



generally not be extracted from the Hipparcos IAD, which span only 3 yrs (see left 
panel of Fig. [191). Similarly, when 07 < 5 mas, as have most barium stars (right panel of 
Fig.[T9l). the Hipparcos IAD are not precise enough to extract the orbital motion. 

When the orbit is not known beforehand, the method makes it even possible to find 
a good estimate for the orbital period, provided that the true period is not an integer 
fraction, or a multiple, of one year. Because parallactic motion has a 1-yr period, any 
orbital motion with a period corresponding to a (sub)multiple of 1 yr will not be easily 
disentangled from the parallactic motion. Therefore, for those periods, the parallax will 
be strongly correlated with the orbital parameters so that the efficiency e (Eq. [251) will 
be small, causing a to be large (see Fig.[20[). 

Inspection of the run of a versus the trial periods P clearly shows that a is minimum 
in the vicinity of the true (spectroscopic) orbital period (see left panel of Fig. [201). The 
period may thus be guessed by looking at the region where a is small. Conversely, 
for stars with no evidence of binarity from a radial-velocity monitoring, neither does 
astrometry find a period range where a becomes small (right panel of Fig. [2T1) . 

Thus, an interesting astrophysical outcome of the algorithm is a list of classical and 
metal-deficient bari um s tars shown to be astrometric binaries by the analysis of the 
Hipparcos IAD lll8Llll4ll . and which were not subject to radial-velocity monitoring so 
far. 




FIGURE 21. Same as Fig. |20]for the P = 2018 d spectroscopic binary HIP 13055 (left panel) and for 
HIP 944 (middle panel) for which P > 800 d may be inferred. For comparison, the right panel displays 
the a statistics for HIP 58596, a star with no evidence for binarity, neither from astrometry nor from 
specti'oscopy l,l&.,20i] . 



6.3. A/i binaries: Comparing Hipparcos and Tycho-2 proper motions 

Widen [fliol] . Wielen et al. lfl2lll and Ka plan and Makarov lHoiil suggested that the 
comparison of Hipparcos and Tycho-2 1 12211123 *1 proper motions offers a way to detect 



Wielen [|T2Cill . Wielen et al. lfl2lll and Ka plan and Makarov lHoiil suggested that the 
comparison of Hipparcos and Tycho-2 1 12211123 *1 proper motions offers a way to detect 
binaries with long periods (typically from 1500 to 30000 d). The Hipparcos proper 
motion, being based on observations spanning only 3 yrs, may be altered by the orbital 
motion, especially for systems with periods in the above range whose orbital motion 
was not recognized by Hipparcos. On the other hand, this effect should average out in 
the Tycho-2 proper motion, which is derived from observations c overi ng a much l onge r 
time span (Fig. [22l) . This method has been used by Wielen et al. UZm, Makarov [|l02|l . 
Pourbaix [|ll3n and Frankowski et al. [llOOn . 
The method evaluates the quantity 



.^obs 



(MhIP - MXyc ) ' W ^ (/IhiP - MTyc ) 



(27) 



where /Ihip and /Ixyc are the vectors of a and 5 components of the Hipparcos and Tycho- 
2 proper motions, respectively, and W is the associated 2x2 variance-covariance matrix. 

Since the above quantity follows a probability distribution function with 2 degrees 
of freedom, it is then possible to compute the probability Prob(;f^ > Xohs)^ giving the 

first kind risk of rejecting the null hypothesis /ixycho = Uuip while it is actually true^ 

To evaluate the efficiency of the method, it has been applied by Frankowski et al. lIlOOll 
to all spectroscopic binary stars from the 5g9 catalogue [jlllll with both an Hipparcos and 
a Tycho-2 entry. Figs. [23] and [24] show that the detection efficiency is very good in the 
period range 1500 - 30000 d for systems with parallaxes in excess of 10 - 20 mas. On the 
other hand, when the proper-motion method detects systems with orbital periods shorter 
than about 400 d, there is a good chance that the system is triple, the proper-motion 
method being sensitive to the long-period pair. Thi s sus picion has been confirmed by 



Frankowski et al. [100], and by Fekel et al. [Ijl, [125 



for the two 'textbook' cases 

HIP 72939 and HIP 88848. Makarov and Kaplan [l 10311 and Frankowski et al. iflOOll 



flagged these two stars as proper-motion binaries, despite the fact that the orbital periods 
known at the time are quite short, 3.55 and 1.81 d respectively. The discovery by Fekel 
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FIGURE 22. The orbital motion (solid line) of a binary, when sampled with a time span short with 
respect to the orbital period, may remain unnoticed, being an almost linear motion (represented by the 
arrows), adding to the proper motion of the centre of mass. On the contrary, when based on positional 
measurements covering a time span much longer than the orbital period, the proper motion will not be 
affected by the orbital motion, which averages out. (From Il20ll ) 



TABLE 2. Two systems whose triple nature was first suspected 
by the proper-motion method and later confirmed by radial- 
velocity studies. 



HIP 72939 
P = 3.55 and 1641 d 
Fekel et al. [125] 





Hipparcos Tycho-2 


Hip. reproc. 


jj.aCos8 (mas/yr) 


-67.5±0.9 -66.2±0.8 


-67.0±1.9 


(mas/yr) 


55 .9 ±0.9 62.1 ±0.9 


63.5±1.5 




HIP 88848 






P= 1.81 and 2092 d 






Fekel et al. [124] 






Hipparcos Tycho-2 


Hip. reproc. 


fiaCosS (mas/yr) 


138.1±1.9 108.5±1.3 


107. 


jUg (mas/yr) 


-18.6±1.7 -25.4±1.2 


-31 



et al. 1 1241112511 that these two systems are in fact triple, with long periods of 1641 and 
2092 d, respectively, came as a nice a posteriori validation of the proper-motion method 
for identifying long-period binaries. The reprocessing of the Hipparcos data, accounting 
for the long-period pair, then yielded values for the Hipparcos proper motions perfectly 
consistent with the Tycho-2 ones, as it should. The situation is summarized in Table[2l 

A very interesting side product of the proper-motion binary detection method is that, 
being totally independent of the star spectrum, it provides an unbiased estimate of the 




FIGURE 23. Eccentricity-period diagram for stars from the 5^9 catalogue with parallaxes larger than 
10 mas. Symbols are as follows: large filled symbols: proper-motion binaries detected at the 0.9999 
confidence level, defined as 1— Prob(;if^ > Xobs) (circles and squares denote those stars which are or are 
not, respectively, flagged as DMSA/G in the Hipparcos Catalogue, i.e., as stars with a non-linear proper 
motion); crosses: systems not flagged as DMSA/G, lying in the confidence-level range 0.99 - 0.99 99; 
open circles: DMSA/G systems not flagged as proper-motion binaries at the 0.9999 level. (From llOOll ) 



binary frequency among the different spectral classes. Table |3]reveals that the frequency 
of proper- motion binaries is cons tant within the error bars (at least in the range F to M). 
According to Frankowski et al. proper-motion binaries represent about 35% of 

the total number of spectroscopic binaries, because of their more restricted period range. 
Thus, the frequencies of proper-motion binaries listed in Table [3] can be multiplied by 
1/0.35 = 2.86 to yield the total fraction of spectroscopic binaries, namely 28.8%. This 
value is close to the spectro scop ic -binary frequency of 30.8 ± 1.7% for cluster K giants 
found by Mermilliod et al. lll26ll . 



6.4. Variability-Induced Movers and Color-induced displacement 

binaries 



The two methods presented in this section are unusual in the sense that they combine 
photometry and astrometry to diagnose a star as a binary. 




FIGURE 24. Upper left panel: Parallax disttibution of Sg9 stars with orbital periods larger than 1500 d. 
Lower left panel: Fraction of those 8^9 stars detected as proper-motion binaries at confidence levels a of 
0.99 and 0.9999. For parallaxes between 10 and 25 mas, about 50% of the S^g stars with periods larger 
than 1500 d are detected, and that fraction becomes larger than 80% above 25 mas. Upper right panel: 
Orbital-period distribution for stars from the S^g catalogue. Lower right panel: Fraction of stars in a given 
bin for which the Hipparcos an d Tyc ho-2 proper motions differ at confidence levels of 0.9999 (dotted line) 
and 0.99 (dashed line). (From iflOOll ) 



TABLE 3. The frequency of proper- 
motion binari es am ong different spectral 
classes, from lIlOOll . is the total number 
of stars of the corresponding spectral type 
in the sample. 



Spectral type 


% 




B 


14.3±13.2 


8664 


A 


15.2±3.0 


15662 


F 


11.4±1.0 


21340 


G 


10.5 ±0.8 


19628 


K 


9.3±0.7 


29349 


M 


9.8±1.3 


4217 


All 


10.1 ±0.4 


103304 



The first category, Variability-Induced Movers, is defined in the Double and Multiple 
System Annex (DMSA) of the Hipparcos Catalogue, where it is known as DMSA/V. 
Stars flagged as DMSA/V are generally large- amplitude variables whose residuals with 
respect to a single-star astrometric solution are correlated with the star apparent bright- 
ness. Such a situation may arise in a binary system, when the position of the photocentre 
of the system varies as a result of the light variability of one of its components: when the 
variable component is at maximum light, it dominates the system light, and the photo- 
centre position is closer to the position of that star [|l27L see also Eq. |2Ql. Pourbaix et al. 
[|l28 1 have shown, however, that many of the long-period variables DMSA/V listed in 
the Hipparcos Catalogue are not binaries, the correlation of the astrometric residuals 



TABLE 4. The 27 stai's keeping their DMSA/V status after the appro- 
priate re-processing by Pourbaix et al. lll28ll . They may be true binaries. 
The column labelled ' Var' lists the variabiUty type (M = Mira; SR = semi- 
regular; L = irregular). 
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with the stellar brightness resulting in fact from the colour variation accompanying the 
light variation of a long-period variable. This colour variation was not accounted for in 
the standard processing applied by the Hipparcos reduction consortia, which adopted a 
constant colour for every star. Thus, the so-called chromaticity effect (for a given po- 
sition on the sky, the position of a star on the Hipparcos focal plane depends on the 
star colour, because of various chromatic aberrations in the Hipparcos optical system) 
was not appropriately corrected in the case of long-period variables. When this effect 
is incorporated in the astrometric processing [by using epoch colours; 'l29], a single- 
star solution appears to appropriately fit the observations for 161 among the 188 stars 
origi nally flagged as DMAS/V in the Hipparcos Catalogue, as shown by Pourbaix et al. 
[fl28i] . The 27 remaining DMSA/V may be true binaries, and are listed in Table IH 

Colour-induced displacement binaries are detected when the position of the photocen- 
tre depends on th e co nsidered photometric band, after correcting for any instrumental 
chromatic effect il30|l . This situation is analogous to the previous one, except that the 
light variability of one component is here replaced by the colour difference between the 
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FIGURE 25. Schematic position of the photocentre in the different SDSS bands. Right panel: For 
double stars, the positions are aligned along the line joining the two stars and are ordered according to the 
central wavelength of the filter The true positions of the stars are represented as a five-branch star symbol. 
Left p anel: For single stars, measurement errors prevent the positions from being perfectly superimposed. 
(From CH) 



two Stars: if the two components of a binary system have very different colours, the re- 
spective contribution of each component to the integrated system light will be different 
in the different photometric bands, and so will be the position of the photocentre. This 
method thus requires accurate astrometry in various filters and is able to detect only bi- 
nary systems with components of very different spectral types. The former condition is 
met by the Sloan Digital Sky Survey rSDSS: ll3lll . which provides positions in the Gunn 
system u,s,r,i, z (F ig. l25l) . The latter condition is met by white-dwarf/red-dwarf pairs. 
Pourbaix et al. [1 132) fou nd 346 pairs among the 4. 1 10^ stars with m, ^ < 21 in the second 
SDSS data releasePBlll with a distance between the u and z positions larger than 0.5 . 
Most (about 90%) of these must indeed correspond to red-dwarf/white-dwarf pairs, as 
indicated on Fig. [261 the remaining 10% being probably early main-sequence/red-giant 
pairs. 



7. OTHER DETECTION METHODS 
7.1. Detecting binaries from rapid rotation 

Another method for finding binaries (among old, late-type stars) involves the identifi- 
cation of fast rotators, since evolved late-type stars are not expected to rotate fast [see e.g. 



133L 11341] . Fast rotation can be ascribed to spin-up processes oper ating either in tidally 
interacting systems [like RS CVn systems among K giants; Il35ll . or in mass-transfer 
systems, through transfer of spin angular momentum. 

In recent years, the link between rapid rotation of old, late-type stars and spin- 
accretion during (wind) mass transfer in binaries has been strengthened, most notably 



2 




u-g 

FIGURE 26. Color-color diagram of the putative binaries (triangles) superimposed on the original 
parent population from the second SDSS data release (contours). The thick lines on the left correspond 
to systems with an M dwarf and three different kinds of WDs. Thin lines are the same, but for a K7V 
star instead of an M dwarf. The short thick solid line in the centre, overlapping the density contours, 
corresponds to A0V/K5III pairs. (From fmll ) 



by the discovery of the family of WIRRing stars [standing for 'Wind-Induced Rapidly 
Rotating'; @]. Originally, this class was defined from a small group of rapidly-rotating, 
magnetically-active K dwarfs with hot WD companions discovered by the ROSA T Wide 
Field Camera and Extreme Ultraviolet Explorer (EUVE) surveys {% Since 



these stars exhibit no short-term radial- velocity variations, it may be concluded that 
any orbital period must be a few months long at least. Moreover, several arguments, 
based on proper motion, WD cooling time scale, and lack of photospheric Li, indicate 
that the rapid rotation of the K dwarf cannot be ascribed to youth. Jeffries and Stevens 
dsl] therefore suggested that the K dwarfs in these wide systems were spun up by the 
accretion of the wind from their companion, when the latter was a mass-losing AGB 
star. The possibility of accreting a substantial amount of spin from the companion's 
wind has been predicted iSlu ,58,. .59.1 by smooth-particle hydrodynamics simulations of 
wind accretion in detached binary systems. A clear signature that mass transfer has been 
operative in the WIRRing syst em 2RE J0357-I-283 is provided by the detection of an 
overabundance of barium [|l36|l . Interestingly enough, the class of WIRRing stars is no 
more restricted to binaries with a late dwarf primary, since in recent years, WIRRing 



systems were found among many different classes of stars: 



A number of barium stars (post-mass-transfer systems involving a K giant and a 
WD ) have been found to rotate fast, within some cases [HD 165141 and 56 Peg; [35, 



13811 clear indications that the mass transfer responsible for the chemical pollution 



of the atmosphere of the giant has occurred recently enough for the magnetic 
braking not to have slowed down the giant appreciably. These systems thus exhibit 
signatures of X-ray activity comparable to RS CVn systems (HD 165141, 56 Peg) 
or host a warm, relatively yo ung W D (HD 165141). HD 77247 is another barium 
star with broad spectral lines llll7ll . It could thus be a fast rotator, but it has no other 
peculiarity, except for its outlying location in the eccentricity - period diagram 
(P = 80.5 d; e = 0.09). In this case, the fast rotation could thus rather be the result 
of tidal synchronization. 

Among binary M (or C) giants, V Hya, HD 190658 and HD 219654 are fast 
rotators (respectively with Vs'mi = 16 km s-K 15.1 ±0.1 km s^^ and 13.6 ± 
0.8 km s^i) [43,|75D. HD 190658 is the M-giant binary with the second shortest 
period known so far (199 d; [p^]), and exhibits as well ellipsoidal variations \J3^. 
Fast rotation is therefore expected in such a close binary with strong tidal inter- 
actions. The radius derived under the assumption of a rotation synchronized with 
the orbital motion is then Rs'mi = 59 Rq, corresponding io R/a\ = 1.16, where 
ai is the semi-major axis of the giant's orbit around the centre of mass of the sys- 
tem il4(^ . The radius deduced from Stefan-Boltzmann law is 62.2 R© [67], which 
implies an orbit seen very close to edge-on (sin/ = 0.948 or / = 71.4°). Adopting 
typical masses of 1.7 M0 for the giant and 1.0 M0 for the companion, one obtains 
R/a = 0.43 and R/Rj^ = 1 ! Thus, the star apparently fills its Roche lobe. 
Binary nuclei of planetary nebulae of the Abell-35 kind (containing only four 
members so far: Abell 35 = BD-22°3467 = LW Hya, LoTr 1, LoTr 5 = HD 1 12313 
= IN Com = 2RE J1255-I-255, and WeBo 1 = PN G 135.6-1-0 1.0) have their optical 



spectra dominated by late-ty pe (G-K) stars 
in the case of WeBo 1 only: Tl4ll [I42L [l43 



with demonstrated barium anomaly 



144|l . but the UV spectra reveal the 



presence of extremely hot (> 10^ K), hence young, WD companions. A related 
system is H P 128 220, consisting of an sdO and a G star in a binary system of 
period 872 d lll45ll . In all these cases, the late-type star is chromospherically active 
and rapidly rotating, and this rapid rotation is likely to result from a recent episode 
of mass transfer. 

Rapid rotation seems to be a comrnon property of the giant component in the few 
known symbiotics of type d' [Tsi 146 1. The evolutionary status of this rare set of 
yellow d' symbiotic systems (SyS), which were all shown to be of solar metallicity, 
has recently been clarified 1II8I1 with the realisation that in these systems, the 
companion is intrinsically hot (because it recently evolved off the AGB), rather 
than being heated by accretion or nuclear burning as in the other brands of SyS. 
Several arguments support this claim: (i) d' SyS host G-type giants whose mass 
loss is not strong enough to heat the companion through accretion and/or nuclear 
burning; (ii) the cool dust observed in d' SyS 11311 is a relic from the mass lost by 
the AGB star; (iii) the optical nebulae observed in d' SyS are most likely genuine 
PN rather than the nebulae associated with the ionized wind of the cool component 



11 14711 ■ d' SyS often appear in PN catalogues. AS 201 for instance actually hosts 
two nebulae fl48]: a large fossil planetary nebula detected by direct imaging, 
and a small nebula formed in the wind of the current cool component; (iv) rapid 
rotation has likely been caused by spin accretion from the former AGB wind like 
in WIRRing systems. The fact that the cool star has not yet been slowed down 
by magnetic braking is another indication that the mass transfer occurred fairly 
recently. Corradi and Schwarz [| 14911 obtained 4000 y for the age of the nebula 
around AS 201, and 40000 y for that around V417 Cen. 
• The blue stragglers were first identified by Sandage lIlSOll in the globular cluster 
M3. Since their discovery, blue stragglers have been found in many other globu- 
lar clusters and in open clusters as we ll. T he frequency of blue stragglers is also 
correlated with position in clusters [see llSlL and references therein]. Depending on 
the cluster, this frequency may either increase or decrease with the distance from 
the centre, thus calling for two different mechanisms for the creation of blue strag- 
glers in clusters. In the central regions, dynamical effe cts in volving collisions may 
enhance the probabilities for creating blue stragglers il52 \. In the outer regions, 
whe re binary systems are most likely to survive, mass transfer is a likely explana- 
tion 1 153 . 154I 1. Despite being initially found in the globular clusters, blue stragglers 
are now fou nd as well in the (halo) field where they are k nown as either 'blue metal- 
poor stars' lll55Lll56lll57ll or 'ultra-Li-deficie nt stars' [ I58,[l591. The very large 
binary frequency among these classes of stars jlSlL ll56Lll57L Il60tl strongly sug- 
gests that field blue stragglers result from mass transfer in a binary, rather than from 
the various merger processes that are currently believed to produce blue stragglers 
in the inner regions of globular clusters. 

Finally comi ng t o the issue at stake in this section, Carney et al. jlSlI) . Preston 
and Sneden (156'] and Ryan et al. [159] find that the field blue straggler stars that 
are binaries have higher rotational velocities than stars of comparable temperature 
but showing no evidence of binarity. Moreover, the orbital periods are too long for 
tidal effects to be important, implying that spin-up during mass transfer when the 
orbital separations and periods were smaller is the cause of the enhanced rotation. 
This conclusion is supported by the abnormally high proportion of blue-metal-poor 
binaries with long periods and small orbital eccentricities, properties that these 
binar ies share with barium stars and related binaries il56ll . Fuhrmann and Bemkopf 
116111 likewise draw a connection between rapid rotation, Li depletion, and mass 
transfer from a companion in rather long -period systems (several hundred days). 



7.2. X-rays 



Several physical processes, which will be discusssed below, are responsible for the 
emission of X-rays in binaries, so that their detection provides a useful diagnostic of 
binarity. In single stars. X-rays are only expected in stars loc ated to the left of the Linsky- 
Haisch line in the Hertzsprung-Russell diagram 1 162l 163 1. separating warm stars with 
hot coronae from cool stars with wind mass loss (Fig. [271). Exceptions to this rule are, 
however, M dwarfs rotating rapidly because they are young, thus producing a dynamo 




FIGURE 27. Comparison of the location in the Hertzsprang-Russell diagram of stars with and without 
X-rays detected by the ROSAT sat ellite (right and left panels, respectively). The 'dividing hne' is clearly 
apparent on the right panel. (From lll63|] ') 



effect which may heat the corona. 

X-rays (corresponding to photons with hv >0.1 keV requiring T > 10^ K) may come 
from: 



• hot stellar coronae heated by a dynamo triggered by fast rotation [|164ll . either 



due to y outh, or to synchronous rotation in a close binary system [as in RS CVn 



systems; HMD- 

• nuclear fusion at the surface of a WD (as in novae or symbiotic stars). In classical 
novae, explosive H-buming (releasing 10^^ — 10^^ erg s^^) occurs, whereas in 
'super soft sources' (as are symbiotic stars), it is probably quiescent H-buming. 

• accretion in binary systems like 

- cataclysmic variables (CVs, consisting of a dwarf star and a WD in a semi- 
detached system with periods of a few hours), which themselves subdivide 
into 

* Classical novae, some can be detected in X-rays outside of eruptions; 

* Dwarf novae, with outburst powered by accretion-disk instability; 

* Polars (AM Her systems), involving highly magnetized WDs around 
which no accretion disk can form because of the very strong magnetic 
field; 

* Intermediate Polars (DQ Her systems), involving highly magnetized 
WDs in a rather wide system which allows some room for an accretion 
disk restricted to the outer region of the Roche lobe; in the inner region, 
no accretion disk can form because of the very strong magnetic field; 

- Algols (which are semi-detached systems consisting of a subgiant component 
and a more massive main sequence component); 

- symbiotic systems, consisting of a giant star and a WD in a system with a 
period exceeding 100 d. However, as we discuss below, the accretion-driven 
nature of the X-rays emitted by SyS has been debated. 



• wind collision, as in SyS. 



7.2.7. X-rays from binaries: The case of SyS 



One of the defining properties of SyS is to host a hot compact star, generally a 
WD, accreting matter from a giant companion [166]. The WD is heated either directly 
by accretion, or indirectly by nuclear burni ng fu eled by accretion [16, and references 
therein]. Symbiotic stars are X-ray sources il67ll . at the level 10-^*^ to 10-^^ erg s^^. The 
physical process emitting these X-rays in symbiotic stars is still debated, though. Direct 
evidence for the presence of an accretion disk, in the form o f continuum fli ckering and 
far UV continuum is not usually found in symbiotic stars [|168L unlike the 

situation prevailing in cataclysmic variables and low-mass X-ray binaries [|l68L Il7l|] . 
Z And is the only SyS where the accretion-driven nature of X-rays makes no doubt, 
because flickering on time scales ranging from minutes to days has been observed 
[|l69ll . Other mechanisms were therefore advocated to account for the X-ray emission 
from SyS, like thermal radiation from the hot component in the case of supersoft X-ray 
sources, or the sh ock f orming in the collisio n regi on between the winds from the hot and 
cool components [|l67l] . Interestingly, Soker [|l72|] even suggests that the X-ray flux from 
SyS exclusively arises from the fast rotation of the cool component spun up by wind 
accretion from the former AGB component (now a WD)! If that hypothesis is correct, 
the properties of X-rays from SyS would be undistinguishable from those of active (RS 
CVn-like) binaries. It is an intriguing hypothesis, because it requires a high incidence 
of fast rotation among SyS, which has so far only been reported for d' symbiotics [see 

above and [iiij. 

The barium system 56 Peg provides an interesting illustration of the difficulty of 
finding the exact origin of t he em itted X-rays. X-rays were detected with the Ei nstein 
satellite by Schindler et al. [173]. In a subsequent study, Dominy and Lambert il74|l 
attribute d the X-rays to an accretion disk around the WD companion. Frankowski and 
Jorissen il38|l re-interpr eted the X-rays as due to RS CVn-like activity, because the 
recent finding by Griffin [OlTSI] that the orbital period of the system is as short as 1 1 1 d 
can only be reconciled with other properties of 56 Peg if it is a fast rotator seen almost 
pole-on! Its X-r ay lu minosity is consistent with the rotation-activity relationship of 
Pallavicini et al. [164], provided that it rotates at a velocity of about 50 km s^^ (Fig.[28l) 
[despite an observed value Vrotsinz of only 4.4 km s~M 134] . thus implying an almost 
pole-on orbit. A low inclination is also required by evolutionary considerations, in order 
to allow the companion of this barium star to be a WD, given the very low mass function 
obtained by Griffin [175]: f{M) = 3.73 10"^ M©. 

Given the small orbital separation, one may still wonder whether there should not 
be as well some cont ribution to the X-ray flux coming from accretion, as advocated by 
Schindler et al. [|l73|] and Dominy and Lambert U2M- The major evidence suggesting 
that there may be a link between mass transfer and X-ray emission in active binaries is 
the correlation between the X-ray luminosity and the Roche-lob e filli ng fraction 72 for 
RS CVn systems noted by Welty and Ramsey lll76ll . Singh et al. 1 177 1 re-examined this 
issue using samples of RS CVn and Algol binaries, and confirm the weak correlation 
found earlier by Welty and Ramsey [176]. However, Singh et al. "regard this as a rather 
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FIGURE 28. X-ray luminosities vs. projected rotational velocities for 56 Peg (encircled crosses, for 
y sin / and for two values of V: assuming orbital corotation or fast rotation) and stars of various spectral 
types and luminosity classes detected by the Einstein Observatory. Different symbols indicate different 
luminosity classes, as indicated in the figure. (Adapted from 1,138,1 and 1,164.1 ). 



weak argument, because both the X-ray luminosities and the Roche lobe filling fractions 
are themselves correlated with the stellar radii in the sample of RS CVn binaries, and 
thus regard this correlation as merely a by-product of the inherent size dependence of 
these quantities." Actually, despite the short orbital period, the filling factor of 56 Peg A 
is not that close to unityO 



7.3. Composite spectra and magnitudes 

SyS are a specific example of a more general class of binaries having composite spec- 
tra, where the spectral features of the two components are intermingled. Spectroscopi- 
cally, these systems correspond to SB2 binaries. Ginestet et al. [|180ll and Ginestet and 



^ Adopting / = 5°, Mi =3 M©, M2 = 0.96 M©, Griffin's orbital elements yield an orbital separation 
A = 152 R0 and a Roche lobe radius around the giant component of 73 Rq, well in excess of the 40 Rq 
representing the stellar radius itself, and translating into 72 = 0.55. 
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FIGURE 29. Left panel: Theoretical zero-age single star and equal-mass binary main sequences for a 
metallicity of Z — 0.02 covering a range of absolute magnitudes My- Also plotted for a range of primary 
masses. Mi, are binary points with M2 = qM\ where q ranges from 1.0 ^ 0.0 in increments of 0.1. 
The point at <7 = 0.5 is an open square. (From llTSll '). Right panel: The (y,Z? — V) color-magnitude 
diagram for low-mass Praesepe members. Symbols are as follows: filled squares: optically -resolved binary 
systems; filled dots: known spectroscopic binaries; filled triangles: triple systems; open circles: 's ingle ' 
stars (although some of these are resolved in the infrared); crosses : suspected non members. (From 11721) 



Carquillat [|181|] provide lists of more than 100 systems with true composite spectra of 
various spectral types and luminosity classes. 

In color-color and color-magnitude diagrams, these stars occupy outlying positions 
and may therefore be easily identified. In colour-magnitude diagrams of clusters for 
instance, binar ies widen the main sequence. This effect has been quantified by Hurley 
and Tout [| 17811 . as shown on the left panel of Fig. [291 A beautiful illustration of this 
effect is provided by the color-magnitude diagram of the Praesepe open cluster where 
the binary sequence is clearly apparent (Right panel of Fig. [29l) . 



7.4. Detecting binaries from the partial absence of maser emission 

Maser emission is observed in many LPVs, and corresponds to stimulated molecular 
emission from an upper level which is more populated than the lower level (population 
inversion). For a review of physical processes related to astronomical masers, we refer to 
the textbook by Elitzur [182]. The maser emission originates from gas in the circumstel- 
lar nebula (fed, e.g., by the strong stellar wind of late-type giants) and involves various 
molecules: SiO, H2O and OH in O-rich environments, and mostly HCN in C-rich en- 




D ID 



FIGURE 30. The maser emission around the supergiant star VX Sgr The solid and dashed lines 
correspond to contours of iso-intensity for the 1612 MHz OH maser line, for gasline-of-sight velocities 
of 6.5 and 1 1.9 km s^^, respectively. The stellar velocity is 5.5 km s^^ (From 1 184ll ) 



vironments. Depending on the energy required to populate the upper level, the maser 
emission will originate from regions closer or farther away from the central star. As dis- 
cussed by Elitzur [182] and Lewis [183], the SiO maser requires the highest excitation 
energies (with the involved levels lying at energies between k x 1730 K and k x 5500 K, 
where k is the Boltzmann constant) and thus originates at the closest distances from the 
star. The H2O maser requires intermediate excitation temperatures (the involved levels 
correspond to E/k between 300 and 1900 K), and operates at intermediate distances 
from the star, whereas the OH maser (subdivided into the so-called main OH maser at 
frequencies 1665 and 1667 MHz, and the satellite line 1612 MHz OH maser) involves 
levels with low excitation temperatures (~ 200 K). Therefore, when the three 0-bearing 
masers operate concurrently in a circumstellar shell, they do so at different distances 
from the central star. A textbook example is provided by the supergiant star VX Sgr 
[Fig. 1301 [l84l]. 

The ro le of binarity in suppressing the maser activity was suspected on theoretical 
grounds [|l85n . since the presence of a companion periodically disturbs the layer where 
maser activity should develop. The systematic absence of SiO masers in symbiotic stars 
except in the very wide system R Aqr jisd Il87|] is evidence supporting that claim. 



On the other hand, very wide systems like o Ceti, R Aqr and X Oph lack the OH 
maser, which involves layers several 10^ R0 away from the star, in contrast to the 
SiO maser, which forms close to the photosphere. There should actually be a critical 
orbital separation for every kind ofmaser below which the companion sweeps through 
the corresponding masing layer 1II86I1 . thus preventing its development (A < 10 AU: no 
masing activity; 10<A(AU) < 50: SiO and H2O masers possible but no OH maser; 
A > 50 AU: all masers can operate). 

A survey of IRAS sources by Lewis et al. [|l88l] reveals that the region occupied by 
OH/IR sources in the IRAS color-color diagram also contains many stars with no OH 
masing activity. Furthermore, detected and undetected OH maser sources have identical 
galactic latitude distributions. These two facts find a natural explanation if binarity is the 
distinctive property between detections and nondetections. 

Lewis H 183(1 [see also 190] has presented a chronological sequence of maser states, 
with increasing mass loss and thus shell thickness. The SiO maser, H2O maser, the 
main OH maser and finally the 1612 MHz OH maser are added one by one as the shell 
grows in thickness. This sequence is supported by the fact that SiO and H2O masers are 
associated with optically identified stars, while the majority of 1612 MHz OH maser 
sources lack an optical counterpart (OH/IR sources). The masers then disappear in the 
reverse order as the shell continues to thicken, with the exception of the 1612 MHz 
OH maser, which remains strong. The SiO maser disappears once steady mass loss 
stops. Thereafter, the main OH maser eventually reappears and grows rapidly in strength 
relative to the 1612 MHz maser in the detached fossil shell, until its ongoing dilation 
causes all of the OH masers to fade away. 

Anomalies in this chronological sequence may be another indication that the star 
has a companion which sweeps through the layer where the lacking maser should 
operate. Therefore , a s earch for binaries among late-type giants involves looking in 
maser catalogues [jl9ll] for 0-rich stars with SiO and OH m asers but no H2O maser, 
since this combination is never encountered in Lewis 11831] chronological sequence. 
RU Her is an example of such an anomaly; it is therefore quite significant that this 
star is suspected to have a composite spectrum [192] and appears in Table |4] as a true 
'Variability -Induced Mover', which are two further indications of its binary nature. 



8. WHAT MAY BE KNOWN ABOUT STELLAR MASSES? 

Binary stars are often considered as the golden way to derive stellar masses. Table [5] 
shows, however, that the different binary-detection methods actually offer only partial 
knowledge of the masses. It is only when a given system appears as both eclipsing and 
SB2 that the individual masses may be derived. We refer to Andersen plQ for a detailed 
review on this topic. 



TABLE 5. A summary of the different properties that may be derived for the various 
classes of binaries. The symbols Ra and d denote the radius of the eclipsed component 
A and the distance of the system from the Sun, respectively. 



Binaries 


Visual 


Astrometric 

/\r>i /\r>Z 


Spectroscopic 


Eclipsing 


MaxP- 




100 yr 


30 yr 


^ 1 yr 


MinP- 


1 yr 


1 yr 


1 d 


1 h 


J? 


yes 


yes yes 


no no 


yes 


a? 


oa + flfi 


oa.b 


a A sin i a^gsin/ 


a 


(units) 


(") 


(") (") 


(km) (km) 


(Ra) 


Masses 


{MA+MB)d-^ 


— Ma.b 


f{M) MA.Bsirr'i 


{Ma+Mb)Ra^ 



AB2 systems (Sirius, a Cen) are rare [189] 

Eel + SB2 (Ex; HD 35155); individual masses may be derived 

9. A LIST OF RESOURCES ON BINARY STARS 

9.1. Combined data 

The best entry point of any search about binary stars is probably the Besan9on 
database http://bdb.obs-besancon.fr/, which provides links to most of the cata- 
logues and databases listed in this section. 

A yearly bibliography about double stars is available at 

|http : 77ad . usno . navy . mll/wds/ dsl/bib' nnnn' . html| 
(where 'nnnn' = 2006 for instance). 

9.2. lAU commissions of interest 

• lAU commission 8 : Astrometry 

[http : / / www . ast . cam . ac . uk/ iau_ comm8/| 

• lAU commission 26 : Binary and Multiple Stars 

|http : / / ad . us no . navy . mil/wds/ dsl . html#iau| 

• lAU commission 27 : Variable Stars 

[EEtp: //www. konkoly.hu/IAUC2 7/ 

• lAU commission 30 : Radial velocities 

[http : / /www . ctio . noao . edu/ science/ iauc30/iauc30 . html] 

• lAU commission 42 : Close Binary Stars 

|http : //www, konkoly. hu/IAUC42/c42"info .html] 



9.3. Visual and interferometric binaries 



Library of links: CHARA (integrates many sources: interferometric, speckle or 
occultation data, WDS, Hipparcos DMSA) 

|http : / / ad . usno . navy . mil/wds/ dsl . html| 

Instantaneous positions: 

- Catalogue of Interferometric Measurements of Binary Stars 

|http : / /ad. usno . navy . mil/wds /int 4 . html] 

- CCDM 

http : / / webviz . u-strasbg . f r/ viz-bin/Vizier ?-source=I /2 7 4| 

- Washington Double Stars (WDS) 

[http : / / ad. usno . navy .mil/wds/wds . html| 

- SIDONIE 

|http : / / sidonie . obs-nice . f r| 

Orbits: 

Sixth Catalogue of Orbits of Visual Binary Stars (Hartkopf & Mason) 

[http : / / ad . usno . navy . mil/wds/ orb6 . html| 



9.4. Spectroscopic binaries 



Spectral families 

- Barium stars llll9ll : 

„http://c dsads. u-strasbg.fr/full/1983ApJS. . .52. .169L| 

- S stars [,1931) : 



[http : //vizier . u-strasbg. f r/viz-bin/VizieR?-source=III/168| 



- C stars 



"http : / /vizier . u-strasbg . f r/viz 


-bin/VizieR?-source=III/22 7| 


Metal-deficient barium and CH stars [ 


195 




T96tl: 



|http: //cdsads .u-strasbg. fr/abs/1990VilOB. .85. 



Is] 



Cataclysmic binaries and X-ray binaries [|197|1 : 

|http : /Tphysics . open . ac . uk/RKcat / (update RKcatV.S, 2007) 

Chromospherically active binaries ['198]: 

http : / / vizi er . u-strasbg . f r/ viz-bin/VizieR?-source=V/7 6| 

Symbiotic stars 

* Ihttp : //vizier .u- str asbg. f r/v iz- bin/VizieR?-s ource =J/A+AS/14 6/4 07[ 

m 



* See as w ell th e series of papers by Carquillat [[199[] . R.E. Griffin [[2001] . 

Ginestet [fisoL fTgii 
Individual radial- velocity measurements 

Ihttp : //www . casleo . gov . ar/catalogo/catalogo . htm] 



Orbits: 

The 9th Catalogue of Spectroscopic Binary Orbits 111 1111 

|http : / /sb9 ■ astro . ulb ■ ac . be/| 



9.5. Astrometric binaries 

Common proper-motion pairs: 

- Luyten LDS catalogue 

[http : / /vizier . u-st Yasbg . f r/viz-bin/VizieR?-source=I/130| 

- Halbwachs 1986 [i20i 

[http : / /vizi er . u-strasbg . f r/viz-bin/VizieR?-source=I/12T| 



- Greaves 2004 1I202I1 



[http': / /vizier . u-strasbg. f r/viz-bin/VizieR?-source=J/MNRAS/355/5"85| 

AjU binaries: 



- DMUBIN database U2M 



|ht tp : / / www . ari . uni-heidelberg . de/ datenbanken/ dmubin| 



- Makarov & Kaplan 2005 ifTosIl 



- Frankowski et al. 2007 ifTooll 
Hipparcos: DMSA/0 binaries 



9.6. Photometric variables 

General Catalogue of Photometric Data 

http : / /obswww . unige . ch / gcpd/ gcpd . html| 

Hipparcos photometry annex 

|http : /7www. rssd. esa . int/SA/HIPPARCOS/apps/PlotCurve . html| 

All Sky Automated Survey 

http : / / archive . princeton . edu^asasT] 
Northern Sky Variability Survey 

|http : //skydot . lanl . gov/| 

OGLE 

http : / /bulge . princeton . edu/~ogle/| 
MACHO 

http : / / wwwmacho . mcmaster . ca/| 

DENIS 

|http : / / www-denis . iap . f r| 

Eclipsing binaries: 

— http : / / cdsweb . u-s trasbg . f r/ cgi-bin/myqcat3 ?eclipsing+binaries| 

- Krakow Observatory 

[http : / / www . oa ■ u j . edu] . 



- Contact (W UMa-type) binaries at Konkoly Observatory 

http : / / www . konkoly. hu /staff/ csizma dia/wuma . html] 

- A catalogue of eclipsing variables [.203,1 
. AAVSO 

|http : / / www . aavso . org| 

• also look for the sites of the French, British, Swiss, New Zealander variable- star- 
observer associations... 
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